Magnetization Oscillation of a Spinor Condensate Induced by Magnetic Field Gradient 
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We study the spin mixing dynamics of ultracold spin-1 atoms in a weak non-uniform magnetic field 
with field gradient G, which can fiip the spin from +1 to —1 so that the magnetization m = p+ — p- 
is not any more a constant. The dynamics of mp = Zeeman component po, as well as the 
system magnetization m, are illustrated for both ferromagnetic and polar interaction cases in the 
mean-field theory. We find that the dynamics of system magnetization can be tuned between the 
Josephson-like oscillation similar to the case of double well, and the interesting self-trapping regimes, 
i.e. the spin mixing dynamics sustains a spontaneous magnetization. Meanwhile the dynamics of 
po may be sufficiently suppressed for initially imbalanced number distribution in the case of polar 
interaction. A "beat-frequency" oscillation of the magnetization emerges in the case of balanced 
initial distribution for polar interaction, which vanishes for ferromagnetic interaction. 
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I. INTRODUCTION 

Since the successful realization of ^■^Na condensate in 
the optical trap with the spin degrees of freedom 
liberated, the coherent spin-mixing dynamics inside a 
spin-1 BEC has been studied intensively [2(-i8|]. The spin 
mixing interaction allows for exchanging atoms among 
spin components but conserving the total angular mo- 
mentum. Two atoms in Zeeman state |0) can coher- 
ently scatter into the states |1) and | — 1) , and vice versa: 
2 |0) ^ |1) -I- | — 1). As one of the most active topics 
in quantum gases, such spin-exchang dynamics was first 
studied by Law et. al. Q and has been observed in 
the way of population oscillations of the Zeeman states 
inside ^^Rb condensates where atoms interact ferro- 
magnetically. A temporal modulation of spin exchange 
interaction, which is tunable with optical Feshbach reso- 
nance, was recently proposed to localize the spin mixing 
dynamics in ^''Rb condensate 

The properties of a three-component {F — 1) spinor 
condensate are first studied by Ho and Ohmi p^ . 
For a spin-1 system, atom- atom interaction takes the 
form V'(r) = (5(r)(co + C2F1 • F2), where r is the dis- 
tance vector between two atoms, and co,C2 denote spin- 
independent and spin-exchange interaction respectively. 
Many predictions are verified experimentally [lll[i2|; and 
the most fundamental property concerns the existence of 
two different phases determined by C2: the so-called po- 
lar (c2 > 0) and ferromagnetic (c2 < 0) states, corre- 
sponding to the F = \ state of ^'^Na and *^Rb atomic 
condensates respectively. The fragmented condensate in 
a uniform magnetic field can be turned into a single con- 
densate state by a field gradient [H, [l3|. More exotic 
ground state phases both in the mean field level and the 
fully quantum many body theory have been extended 
to condensates with hig her spins |15l - [l7| and recently to 
spinor mixtures p^l - [20 |. 



In this paper we study the dependence of the spin dy- 
namics on a small magnetic field gradient, which prac- 
tically provides a process to flip the spin between -fl 
and —1 states thus turns a fragmented condensate into 
a coherent one. We adopt the mean-field approximation, 
in which a spinor condensate is described by a multi- 
component vector field. It has provided surprisingly good 
descriptions for most properties of the spinor condensate 
as evidenced by the experimental verification of many 
predictions . As the spin flipping term is considered in 
the condensate, the system magnetization exhibits obvi- 
ous macroscopic oscillation similar to the Josephson os- 
cillation of a scalar condensate in a double well [2l| and 
we find that the dynamics of spin-0 component po may be 
greatly suppressed in the case of polar interaction. This 
provides us an intriguing tool to manipulate the atomic 
population in spinor condensate. 



II. THE EFFECTIVE HAMILTONIAN FOR THE 
SYSTEM 

The many-body Hamiltonian of N spin-1 atoms of 
mass M in a uniform magnetic field reads 



H 
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E( 

fe=i 



2M 



- Vtrap + 7B0 • Ffc) + V V"^„t (rfe - r;). (1) 



k<l 



The atoms are loaded into an optical trap Vtrap and Vint 
denotes the collisional interaction between atoms. 
and Ffc are the momentum operators and spin-1 opera- 
tors of the fc-th atom, respectively. We consider a field 
gradient G as has already been applied in the MIT ex- 
periment [llj |. and to be more specific we just replace Bq 
in ^ with 



B(r) ^ BoTb ^ Bo[z. + G(a;x - zz)]. 



(2) 
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We choose the local field direction re as the spin quan- 
tization axis [13I I by performing a unitary transforma- 

N . 

tion J7 = n e^i"^'^''^'^'' on the Hamiltonian ([T]) where 

k=l 
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n = z X re =Gxy. The contact interaction Vi„t is spin 
conserving, hence it is invariant under the transforma- 
tion. On the other hand, an additional Berry phase as- 
sociated with the local change of basis emerges in the 
momentum term. Take the single particle hamiltonian as 
an example, h = ^ + Vtrap + l'B-F, and U — e^T^'^^^^, 
we find that 



= -2GFyP^ + G'^F^ 



(3) 



where the operator A — GFy. The Zeeman energy term 
is transformed into 221 



U^B-FU = BoF^[Gxcos{Gx)-{l-Gz)sm{Gx) 
+ BqF, [Gx sin (Gx) + (1 - Gz) cos (Gx) 



(4) 



For small field gradient G, we approximate the trigono- 
metric functions up to order of O(G^) in accordance with 
Eq. ©, cos(Gx) ~ 1 + G'^x^/2,sm{Gx) ~ Gx, which 
leads us to 

U^'B ■ FU ^ Bo[F, + CG + DG^ + 0{G^)] (5) 

where the operator G = —zF^, D = xzF^ + x'^Fz/2. 

The many-body Hamiltonian is finally transformed 
into 



H, 



eff 



k=l 



+Vtrap + 7^0-^^ + T^oGG^' + iB^G^D^} 

k<l 

The atomic interaction takes the form [§, [loj 

V^»„t(r) = (co-}-C2Fi-F2)J(r) 

and Co = '^'Kh^{aa + 2a2)/3M, C2 = AT:h^{a2 - ao)/3M. 
The representation of the Hamiltonian in the second 
quantized form reads 



+ Vtrap + Po{Fz)^J + PqGG.j + PQG'^D,j]i'j 

-h|*1*I,F., -F.,,,*,,*,} 



(7) 



where repeated indices are to be summed over and 
*(r)(*t(r)) is the field operator that annihilates (cre- 
ates) an atom in the i-th hyperfine states with i ~ 1, 0, —1 
at location r, and po = 7^o with 7 the gyromagnetic ra- 
tio of the bosonic atoms. 



The dynamics of the condensate components reveal a 
rich coupling between the spin and spatial degrees of free- 
dom resulting in a variety of interesting phenomena, in- 
cluding spin mixing, spin domain formation and spin tex- 
tures [4]. The internal and external dynamics are both 
very sensitive to the external magnetic fields and field 
gradients and they can be decoupled under certain con- 
ditions, in particular, when the available spin dependent 
interaction is insufficient to create spatial spin structures 
in the condensates. This occurs when the spin healing 
length is larger than the size of the condensate, which al- 
lows us to focus on the coherent spin mixing oscillation of 
the spin populations. Practically we choose a proper field 
gradient G to induce an energy in the same magnitude of 
spin interaction term C2. The gradient will fiip the spin of 
atoms in the condensate but keep the three components 
still miscible and free of spin texture (G < 2cm^^ as in 
the MIT experiment As a result we can still safely 

adopt the single spatial mode approximation (SMA) in 
the following |M,13,|2l- 

We take 



*1 = ViV^(r)a| 



(8) 



with <j>{r) is defined by the Gross-Pitaevskii equation 
through the spin-independent part Hq 
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(j)^Hacp^ lic^ (9) 



where N is the total number of the atoms and /i is the 
mean field energy or the chemical potential. The a|(ai) 
is the spin component operator that annihilates (cre- 
ates) an atom with spin i (i = 1,0, —1) . Substitute the 
ansatz ([8]) into the Hamiltonian ([7]) and neglect the spin- 
independent part, our model then reads 



Heff = -e(a|a_i+aLiai)+ea5ao+C2F2 - pF^ (10) 



with 



F2 

F, 



F^ + {F+F^+F^F+)/2 
V2(a|tiao + aja^i) 



a\ ai 



(11) 



The parameter e — G'^/AM characterizes the spin- 
fiipping process induced by the field gradient, and 
the spin interaction parameter is scaled as c'2 — 
(C2/2) /(ir|$(r)| and we keep the original notation C2 
for ease of representation. The term Cij vanishes due to 
the fact that in ground state is a symmetric function 
(/ dv(f)*x4> = and / dr(j)* z4> = 0, etc.), so does the term 
Px{Fy)ij because 

j dr(P*P,(b =^J d-^^*[^^ W = 0. 

The term Dij amounts to a shift of the linear Zeeman 
energy p = po + p with p = pqG'^ J dr(j>*x'^(j). 
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The e term in Hamiltonian ([10)) denotes the process 
that flips the spin from 1 to —1 or vise versa, and it 
plays the same role as the hopping term of scalar conden- 
sates in a double well. We want to emphasize that this 
term in the form of (x\{Fy)ijaj induces the oscillation of 
the z-component magnetization m and goes against the 
quadratic Zeeman effect term al{F^)ijdj in Refs. [3,[3|> 
which instead adheres to the constancy of the magneti- 
zation m and keep the system in the polar phase 11]. In 
a uniform field, the z-component magnetization m is a 
constant, and the system can be described by two canon- 
ical conjugate variables: the population on spin-0 com- 
ponent po and the relative phase = 6i + — 29o. The 
system can be finally reduced to a nonrigid pendulum 
model 4]. Our model system conserves only the total 
number of atom and another relative phase 9' — 9i — 9^i 
arises. We will mainly study the dynamics of the pop- 
ulation on spin-0 component pait) — NQ(t)/N, and the 
z-component magnetization m(t) = {Ni{t) — N-i{t))/N. 
The populations on ± components are related to po and 
TO through pi + p^i + pq ~ 1 and m — pi ~ p^i. 



III. THE SEMICLASSICAL MODEL 

Instead of considering a Hilbert space of the Fock 
states |A''i, TVq, iV_i), where the population dynamic can 
be describe as 

po^ {iJ^\e"'"'alaoe''f'"^\4>^) /N 

with an initial state which can be taken as one of the 
basis in the Hilbert space, here we consider the conden- 
sate to be in a coherent state, associated with a macro- 
scopic wave function with both magnitude and phase. In 
the study of spin mixing dynamics 3, 7, 8], it is custom- 
ary to replace the operator a\{ai) by c numbers a*(ai), 
and the coherent state is analogous to a classical field 
of complex amplitude with a definite phase in each spin 
component 

1$) = |ai,ao,a-i> - \^ie''\ ^ae''\ ^/N~ie''-') 

Semiclassical equations of motion can be derived from 
Hamiltonian (|10p as 

ihhi — 2c2{Fzai + F^aa/Vi) — ea^i — pai 
ihao = V2c2(#+ai -t-F_a_i) -feao (12) 
ihci-i — 2c2(— i^zC-i + /+ao/\/2) — ea+i + pa_i 

where the quantities F±. , F^ are c number counterparts 
of the operators F± and F^ respectively. The neglected 
spin-independent part of the Hamiltonian ([7]) give each 
of the three equations in Eq. (|12l) a constant energy 
shift p that can be trivially eliminated by changing ai 
to Oje"*''*/''. For simplicity we rescale the time as t 
\ci\tlh.. The three coupled Gross-Pitaevskii equations 
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FIG. 1: (Color online) The dependence of the dynamics of 
poit) (red sohd) and m(t) (blue dashed) on the parameters 
of e at fixed values of p = 0, C2 = 1 (in units of |c2|), and 
e = (a), 1.45 (b), 1.49, 1.50 (c), 1.51 (d), 1.55 (e) and 2.25 (f). 
Fig. 1(c) shows the critical transition parameters of po{t), m{t) 
with e =1.50 (black dashed dot line). Time is in units of 

\c2\t/h. 



(|T2)l for the interacting condensate amplitudes a* {at) 
describe the dynamics in terms of the inter-component 
phase difference and population imbalance. 



IV. RESULTS AND DISCUSSION 

First, we consider an initial distribution with imbal- 
ance between |-fl) and | — 1) components, i.e. |$(0)) = 

0, ^iV/2e*^o , ^iV/2e*^-i^ and illustrate the dynamics of 

population po{t) and magnetization m(t) for both C2 > 
and C2 < 0. In order to give prominence to the effect of 
the e term, we consider p — first, and set the phases 
9q = 9^1 = 0. This initial imbalance iVi(O) - iV_i(0) = 
—N/2 provides a "Junction voltage" and the magnetiza- 
tion oscillation was induced by e. In Fig. [l] we show 
the solutions of equations for C2 > and illustrative 
parameters e = 0, 1.45, 1.49, 1.50, 1.51, 1.55 and 2.25, in 
the unit of |c2|, respectively. We find that at the very 
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FIG. 2: (Color online) The dependence of the dynamics of 
pQ{t) (red solid) and m{t) (blue dashed) on the parameters 
of e at fixed values of p = 0, C2 = —1 (in units of |c2|), and 
e = (a), 0.45 (b), 0.49, 0.50 (c),0.51 (d), 0.85 (e), and 1.50 
(f). Fig. 2 (c) shows the critical transition parameters of m{t) 
with e =0.50 (black dashed dot line). Time is in units of 

\c2\t/h. 



beginning when e is small the magnetization oscillates 
with small amplitude around an equilibrium above the 
initial value of m(0) = —0.5 as in Fig.([T|3), which is 
analogous to the "macroscopi c q uantum self-trapping" 
effect in double well system [2l|. Meanwhile the dy- 
namics of po(t) experiences a crossover from sinusoidal 
to non-sinusoidal oscillation, with the population pait) 
averaged over time changing from less than the initial 
value po(0) — 0.5 to larger than it. As e increases, there 
is a critical transition for e = 1.50, black dashed line in 
Fig-iHh), then the oscillation extends to the range be- 
tween —0.5 and 0.5. Accompanied by the arising of the 
" Josephson tunneling" [2l| of the magnetization, the dy- 
namics of pq has been sufficiently suppressed in Fig.(ITf). 
The coherent scattering of the internal Zeeman compo- 
nents 2|0) ^ |1) + | — 1) was suppressed by the e term 
with the process |1) ^ | — 1) • The critical behavior de- 
pends on e, as can be easily found from the energy con- 
servation and the extreme point for the minimization 
of the energy. Considering an arbitrary wave function 
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FIG. 3: (Color online) The dependence of the dynamics of 
po(i) (red solid) and m{t) (blue dashed) on the parameters 
of p at fixed values of e = 1.55, C2 = 1 (in units of |c2|), and 
p = (a), 0.2 (b),0.4 (c), 0.6 (d). Time is in units of \c2\t/h. 
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FIG. 4: (Color online) The dependence of the dynamics of 
poit) (red solid) and m{t) (blue dashed) on the parameters of 
p at fixed values of e = 0.55, C2 = — 1 (in units of |c2|), and 
p = (a), 0.2 (b),0.4 (c),0.6 (d). Time is in units of \c2\t/h. 



1$) = Ixe*^! , ye*^", 2;e*^-i) , the relative average energy 
of the system when p = can be described as 



E ^ C2[2y^{x^ + z-" + 2x z cose) + {x^ ~ z^y] 



~2exz cos 6' 



(13) 



withe* = e++e^i-2eo, and e' 



For C2 > 0, the 



critical point favors that y = y/2/3, x ^ z — -^/l/G, 6 — 
TT, 0' = TT. According to the energy conservation condition 
Ec — Einitiai, we can get the critical value ec/c2 = 1.50. 
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FIG. 5: (Color online) The dependence of the dy- 
namics of po{t) (red solid) and m{t) (blue dashed) 
on the phases of the initial condition |<1?(0)) = 

|y]v73e'*i,y]v73e^*o,yiV73e'*-i^ with parameters of 

C2 = 1, e = 0.55, p = 0.6 and Oi ^ do ^ O-i = for (a), 
C2 = 1, e = 0.55, p = 0.6 and Si = O-i = 7r/2, Go ^ tt for (b), 
C2 = -1, e = 0.55, p = 0.6 and Oi = 6l_i = So = for (c), 
C2 = -1, e = 0.55, p = 0.6 and 6*1 = 6l_i = 7r/2, So = tv for 
(d). Time is in units of |c2| t/ft. 



effect of the parameter p, while we choose a fixed value 
of e = 1.55 (in the unit of |c2|) for C2 > and e =0.55 for 
C2 < as shown in Fig.([3|) and Fig.dl]). 

We find that the p term acts as a switch for the dynam- 
ics of the magnetization m{t), and the cases for C2 > 
and C2 < are quite different. For C2 > in Fig.([3|), 
the fast oscillation of m{t) is modulated by a beat fre- 
quency, and, when the amplitude of the envelope function 
reaches its maximum the population of po is completely 
suppressed to zero. For the C2 < case, the oscillation 
of 7n{t) is also induced by p, but the magnetization is 
always positive. The dynamics of po shows completely 
anharmonic behavior with the amplitude first enhanced 
then reduced as p increases, and no beat frequency mod- 
ulation of the m{t) occurs. 

However, specific features of quantum nature cannot 
be addressed satisfactorily within a mean- field treatment. 
In an early experiment [5| on an F = 1 ^^Rb condensate, 
atoms all prepared initially in the state |0, N, 0) are ob- 
served to exhibit a damped oscillation accompanied by 
large fluctuations during the spin-mixing evolution. The 
dynamics of a polar initial state |0,A^, 0) is trivial, i.e. 
the population po remains a constant, within a mean-field 
treatment, but many-body quantum dynamics shows in- 
teresting damped oscillation . The presence of a field 
gradient will enlarge the Hilbert space in quantum treat- 
ment due to the failure of the conservation of magneti- 
zation m and the related calculation on this feature will 
be published elsewhere. 

Finally, let's consider the effect of the phase difference 
between the three components. Fig. ([5]) shows the dynam- 
ics in the presence of an initial phase difference for both 
C2 > and C2 < cases. We find that the influence of 
the phase difference for both cases are obvious. The beat 
frequency oscillation of the magnetization remains in the 
case of C2 > with a shift of the envelop center, but for 
the C2 < the phase difference changes the amplitude of 
m{t) which extends down to the negative part of the axis 
and the dynamics of po is also altered drastically. 



For C2 < 0, we illustrate the dynamics for different pa- 
rameters in Fig.dH) with e = 0,0.45,0.49,0.50,0.51,0.85 
and 1.50, in the unit of |c2| respectively. We flnd that 
the oscillation of m{t) here is almost the same as in the 
polar interaction case C2 > 0. On the other hand, po{t) 
is not suppressed at the beginning (Fig.l^b)), instead, it 
was enhanced. When e > 0.5, and the oscillation of po is 
still active until e reaches 1 or even large value. This fer- 
romagnetic feature is quite different from the C2 > case. 
The latter case shows some "repulsive" effect between the 
Zeeman components. The critical value tc/ {—C2) — 0.5 is 
derived analytically through the equation (IT3)) with the 
extreme point y = Q,x = z = a/1/2, 9 = 6' — Q. 

Next, we consider a balanced initial distribution with 

|$(0)) = vW3e**\ \/W3e*^°, \/W3e'^"') and 6^ = 
^-1 = ^^0 = 0, where the "Junction voltage" between 
the |1) and | — 1) components vanishes. We consider the 



V. CONCLUSION 

The dynamics of spin-1 BEG in a nonuniform mag- 
netic field is studied with the emergence of an additional 
spin-fiipping term induced by the field gradient, which 
has an effect to reverse the spin from -1-1 to —1 and 
vise versa. Due to this spin flipping process the system 
magnetization m{t) is not a constant any more, instead, 
it shows characteristic oscillation identical to that of a 
scalar BEG in double well. Meanwhile, the dynamics of 
po was greatly altered. We present the dynamics of pQ{t) 
and m{t) under different initial conditions and the effect 
of phase difference is also shown. We find that the results 
for the polar (c2 > 0) and ferromagnetic (c2 < 0) cases 
are quite different. The small magnetic field gradient is 
chosen properly to give rise to the fiipping of the spin 
between -|-1 and —1 but still keep the three components 
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misciblc. These results highlight the possibility to ma- 
nipulate the coherent dynamics of the spinor condensate 
with a field gradient, which is accessible to the current 
experimental techniques. 
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